Electroweak Baryogenesis from Late Neutrino Masses 
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Electroweak Baryogenesis, given a first order phase transition, does not work in the standard model because 
the quark Yukawa matrices are too hierarchical. On the other hand, the neutrino mass matrix is apparently not 
hierarchical. In models with neutrino mass generation at low scales, the neutrino Yukawa couplings lead to large 
CP-violation in the reflection probability of heavy leptons by the expanding Higgs bubble wall, and can generate 
the observed baryon asymmetry of the universe. The mechanism predicts new vector-like leptons below the TeV 
scale and sizable /j, — > e processes. 
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The neutrino physics has been undergoing a revolutionary 
progress in the past several years. The recent observations 
of neutrino flavor conversion in solar, atmospheric and reac- 
tor neutrino experiments provided firm evidence for physics 
beyond the Standard Model (SM) (see e.g. [[]], §| and refs. 
therein). The new data on the distance/energy-dependence 
of atmospheric neutrinos |||] and the recent spectrum anal- 
ysis from the reactor experiments [Q] indicate oscillatory 
behaviour. This strongly favors the presence of tiny but non- 
zero neutrino masses. 

The most popular explanation for the origin of neutrino 
masses is the seesaw mechanism In its minimal form, 
right-handed neutrinos are introduced with lepton-number 
violating Majorana mass terms, A-^NN, as well as Yukawa 
couplings to the SM lepton doublets I and the Higgs h, 
Y L hNl. The oscillation data requires A L ~ 10 14 GeV if 
Yl ~ 0(1). Moreover, the seesaw mechanism naturally 
provides a way to account for the observed baryon asymmetry 
of the universe (BAU) §|, tib/s ~ 8 X 10~ n , where tlb/s 
is the baryon to entropy ratio. The out-of-equilibrium decay 
of the right-handed neutrinos creates a lepton asymmetry, 
which is partially converted to a baryon asymmetry by the 
electroweak sphaleron process (leptogenesis) 0. Thermal 
leptogenesis typically requires, however, Al ~ 10 9 GeV ^ 
with somewhat small Yukawa couplings Yl and large hierar- 
chies in the right-handed masses. 

While the seesaw mechanism is very appealing theoreti- 
cally, it is unlikely that it will be subject to a direct experimen- 
tal test in the near future. It is important to explore other possi- 
bilities for the origin of neutrino masses. One such example is 
the late neutrino mass framework that induces small neutrino 
masses due to a low scale of symmetry breaking [g]. When 
this symmetry is broken, say by a set of symmetry breaking 
VEVs, ((f)) = /, the neutrinos acquire masses from operators 
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We stress that this does not depend on the details of the 
model, or whether one uses global [^j or gauge [ flOj , lyj [l2|] 
symmetries. The strongest limits on / arise from big bang 
nucleosynthesis plj, |l2| , p"3| , |i4| l and from observation of 
supernova neutrinos [|15||: for n = 1, / > 10 keV, while 
more powerful limits apply to higher n. It is remarkable that 



new physics at such a low scale is not excluded by direct 
experimental data. 

In this letter we show that the late neutrino mass framework 
can naturally realize leptogenesis at the electroweak (EW) 
phase transition [ [lq l . Heavy vector-like leptons that give rise 
to the operators (]l]) bounce off the expanding bubble walls 
with 0(1) Yukawa couplings, and acquire a large asymmetry. 
They quickly decay to the standard-model leptons and the 
sphaleron process partially converts their asymmetry to the 
baryon asymmetry. Thus our mechanism leads to various 
predictions that can be tested by near future experiments. Here 
we focus on the qualitative features of our scenario while 
a more detailed analysis will be presented in a following 
publication. We simply assume that the phase transition is first 
order, and focus on the size of the baryon asymmetry as well 
as phenomenological constraints on the model. We briefly 
comment on the origin of the first-order phase transition 
towards the end of this paper. 

The Model and Mechanism. As well as introducing a 
low scale /, theories of late neutrino masses introduce flavor 
scales that are much lower than the scale A^ of the seesaw 
mechanism. One economical possibility, that we explore 
in this letter, is that these flavor scales are all of order the 
EW scale: M f ,Al ~ v = (h). In this case the non-SM 
states that generate the operators of (jl|) have masses of order 
v and are available to take part in EW baryogenesis. The 
flavor symmetry breaking scale becomes / ~ t^TO^/v) 1 /™ ~ 
(m u , 100 keV, 30MeV...) forn = 1, 2, 3... Then = lease is 
excluded by BBN and supernova constraints. Large n theories 
may be preferred in the sense that the scale / grows and 
requires less protection. In this letter, we describe a simple 
n = 2 theory that illustrates our baryogenesis mechanism. 
We focus on the Majorana case since, as shown below, the 
Dirac case is disfavored by the direct experimental data for 
the parameter range that produce enough baryon asymmetry. 

The minimal model for late Majorana neutrino masses is 

C V = Y hNL + M L C L + y <f>L c l + — NN + h.c. , (2) 

where L, L c are vector-like lepton doublets, and the couplings 
Y, y, and the masses M, Mm are 3x3 flavor matrices. We 
consider all eigenvalues of M to be comparable, the same for 
Mat, and those of Y to be of 0(1). This is suggested by the 
lack of hierarchy (anarchy [f7|) in the neutrino masses. We 



refer to their eigenvalues as M,Mn, and Y . There may be a 
moderate hierarchy in the eigenvalues of y: y^. Note that there 
is a a leptonic analogue of the rephasing invariant Jarlskog 
determinant even for just two generations (Nq — 2) [OX 
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Below the scales M and Mn, neutrino masses are de- 
scribed by the operator (y 2 Y 2 / M 2 MN)l 2 h 2 4> 2 , correspond- 
ing to the n = 2 case of (^), so that the flavor symmetry 
breaking scale is predicted to be 
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A low / ~ 100 keV would imply that the states contribute 
to the energy density of the universe during BBN. However, 
our baryogenesis mechanism will require M, Mn somewhat 
larger than v, and lepton flavor violation will constrain yi to be 
somewhat small. Hence we expect that the <p states are heavier 
than 1 MeV and decay to neutrinos before BBN. 

Let us now describe our main mechanism. With M ~ v, 
during the electroweak phase transition L particles and their 
CP conjugates, L, are reflected differently from the Higgs 
bubble wall. This is due to the presence of unsuppressed 
CP violating phases in Y, M, and Mn (for a related idea 
see JT^ , 19 1). Thus an asymmetry in L is induced in the 
region just in front of the wall. As shown below, the size 
of this asymmetry is expected to be of order Jl/M n . We 
assume the L — > Nh decay process is kinematically forbidden 
(i.e., Mn > M) so that the asymmetry is transferred to 
SM leptons via L — > l<f> decays. The presence of the 
asymmetry in the SM lepton doublets biases the sphaleron rate 
to induce a B production in the vicinity of the wall. When 
the expanding bubble passes over this region, the sphaleron 
processes decouple, freezing in a B asymmetry. Outside the 
bubble the sphaleron rate Ts P ~ a^yT c , T c being the critical 
temperature, is much slower than other dynamical scale near 
the bubble wall [p0[]. Thus the baryon asymmetry could be as 
large as 
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where ~ 100 is the number of relativistic degrees of 
freedom, and anarchical neutrino masses suggest that Jl /M n 
is of order unity. Below we study what other factors might 
suppress the baryon asymmetry. 

Estimating the Baryon Asymmetry. To have a semi- 
quantitative estimation of the resultant BAU, we apply the 
thin wall approximation (the validity of this approximation 
depends on the details of the mechanism which produces the 
1st order phase transition [^o|). We first estimate the density 
difference between L and L, til, induced by the reflection 
asymmetry pit] 
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FIG. 1: Perturbative calculation of the reflection 
coefficients which picks up the Jarlskog invariant Jl = 
QmTr[(M N Yy (Y*Y T M N M* N M N Y)\ from the one-particle 
cuts in the amplitudes. 



where v w ~ 0.1 is the wall velocity and tio(uj) — l/(e"/ Tc + 
1) is the Fermi-Dirac distribution. The difference between 
N and L momenta for a given energy, Ap = p L — pn, 
is large, due to 0(1) mass differences among M and Mn- 
This is welcome because, in the SM, Ap arises only due 
to the electroweak thermal corrections to the masses, and is 
suppressed by aw [pT|]. The reflection asymmetry A(ui) is 

given by A(w) = Tr (r nl Rnl - R nl Rnl) , where 
Rnl is the reflection coefficient for N — > L and the bars 
stand for the CP conjugated process. 

In order to calculate the reflection asymmetry from the 
bubble wall, we have computed the Green function for our 
model in the simplifying limit Mn 3> M. Using a pertur- 
bative expansion in Yv/T c , M/u), and for Nq = 2, we find 

A(uj) - {^t^) 2Ng Jl where l T - 10/T is the mean free 
path for the leptons [ pc| ] (see Fig. |l]). Below we find that 
phenomenological constraints favor somewhat heavier RH 
Majorana masses, Mn ^ Yv. In this case a suppression in 
the reflection asymmetry is expected since the heavy incoming 
particle is hardly affected by the potential barrier. This 
effect cannot be captured using an expansion in Mjy/w. We 
estimated the corresponding 1 /Mn suppression by analysing 
a single generation reflection problem which can be solved 
analytically. We indeed found that the reflection amplitude is 
further suppressed via (Yv / Mn) 2Ng ■ Consequently, in the 
relevant region of parameter space the reflection asymmetry 
is given by 
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Below we use the estimate (^ even in the case of interest 
where Yv/T c and M/u are of order unity. 

The rate for baryon production is approximately dn-Q /dt ~ 
3T c 2 r S pAF [|0|, |2j]. AF is the free energy difference 
between two neighboring zero field-strength configurations, 
for which AB = AL = 3 (provided that T(L l<f>) is not 
much slower than other thermalization rates; see below). AF 
is calculated in the presence of a hypercharge density riy ~ 
— 7j_l/2 [p2j]. (More correctly, one should define a global 
approximately conserved charge orthogonal to hypercharge 
B' = B — xY where in our case x = 1/7 [p2[]. We verified that 
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FIG. 2: n B /s as a function of _M N /T C , for £ w ~ 20/T c , 7V G = 2 
and A(w) = (T c /™) 4 6)(lj - M N ). 



this will hardly modify our results). In our case we find that 
AF ~ ri L /T c 2 . The BAU is obtained by integrating dn B /dt, 
which we estimate via n-Q ~ dn-Q /dt x £ w /v w where £ w is 
the typical penetration length for the non-zero global charge 
which flows through the plasma in the unbroken phase. With 
fast massless leptons one expects, from estimation of energy 
loss in the plasma [g, T c £ w = 0(100) @ g]. In our 
case, with semi-relativistic leptons, the penetration distance is 
shorter even though the energy loss rate is very low, because 
the massive leptons rather quickly lose their directionality 
(away from the wall). For instance, elastic scattering of a 
lepton of mass 4T C with a plasmon that carries a perpendicular 
momenta, T c , roughly results in 25% momentum loss in the 
original direction of motion for the lepton. Thus in our case 
we estimate £ w ~ N co \\£t where iV co u is the average number 
of collisions that a lepton undergoes before its directionality 
is lost. Below, we assume N co \\ = 2. 
Using Eqs. (plpb we obtain 
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Our next step is to identify the dependence of the resul- 
tant BAU on the model fundamental parameters. Assum- 
ing that CP violation is maximal and taking for simplicity 
(Yv) 2 £ T /T c ~ 1 we find A(w) = (T c /w) 4 (9(w - M N ). 
In addition, fixing iV co ii — 2 and £t ~ 10/T C we can 
numerically compute the resultant BAU as a function of 
Mn/T c . As Mm/T c increases, the BAU rapidly decreases 
due to both a Boltzmann suppression factor [see Eq. (|])] 
and also a polynomial one (^. In fig. || we plot Ub/s as a 
function of M/T c . We find that the observed asymmetry can 
be accounted for when 

M N < 4 x T c . (9) 

In the Dirac case Ng = 3 and the resultant asymmetry is 
further suppressed. 

Direct Constraints and Tests. In the following we briefly 
discuss the direct phenomenological constraints on our sce- 
nario and discuss several ways to directly test our model in 
the near future. 

Electroweak precision measurements- Our model requires 
additional fields, namely heavy lepton doublets charged under 



the SM gauge group. Because they are vector-like, the S 
parameter is hardly affected. The Yukawa couplings to the SM 
Higgs, Y, breaks the custodial isospin symmetry and therefore 
modify the T parameter. The additional contribution to the 
T parameter from the three vector-like lepton doublets, Tl, 
is similar to a single extra vector-like top quark [ p4j ] Tl ~ 
where T ( SM ~ 1.2 is the SM-top contribution to 
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T [[!]] and m t is the top mass. Requiring the Tl ^ 0.2 we find 
the following lower bound 
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Furthermore, in the Dirac case, the light RH neutrinos are a 
linear combination of N and L c , and hence the Z can decay 
invisibly to these states. The rate is proportional to the quartic 
power of the corresponding mixing, M4 . Using the 3a range 
(to allow for at least three neutrinos) J2~5| ] we find that 3.01 
neutrinos are allowed. For three extra generations this implies 
M > 4 x Yv . (This constraint is absent in the Majorana case.) 
Lepton flavor violation- The Yukawa couplings to the SM 
leptons, y, are generically not aligned with M, Mm and Y 
and hence induce lepton flavor violation. This will contribute 
to processes such as fi — > e conversion which are highly 
constrained by experimental data ^ |2^]. In our model the 
SM leptons couple to the additional fields only through y, the 
Yukawa coupling to <j>. To avoid these constraints naturally, a 
mild hierarchical structure is required for y. For example, in 
the appropriate basis, y ~ diag (10~ 2 , 10 _1 , 10 _1 ) , which 



is consistent with the neutrino flavor parameters 1 27 1 provided 



that M/v,M, y are anarchical. The smaller the eigenvalues 
of y, yi, however, the smaller is the decay rate L — > lcj>. 
Comparing the decay rate into SM doublets to the thermaliza- 
tion time scale we find T L ( w /v w ~ lOyf M/T c ~ 0.25 - 1 
which implies further suppression in the resultant BAU (this, 
given our crude estimation see fig. ||, is still consistent with the 
observed value). We thus find that there is a tension between 
producing sizable BAU and suppressing the contribution to 
lepton flavor violation. Therefore our mechanism typically 
predicts that the rates for fi — > e processes are within the 
reach of near future experiments. In addition, given the above 
ansatz for y, contributions to neutrinoless double beta decay 
are evaded [[T|] . 

Collider physics- From Eqs. (^ 10), our scenario predicts the 
presence of SU(2) doublets with masses below the Oil TeV) 
range. These contain charged particles, with masses above 
lOOGeV from current bounds [[!]], that may be produced and 
detected at the LHC, or better at the ILC. A clear signal at 
the LHC is expected if the vector-like leptons are sufficiently 
light, via production of L ± L° which decay to a single, 
energetic SM charged lepton. If the Higgs boson is heavier 
than the new particles it can also decay invisibly to and L. 

Note that the presence of new lepton doublets is crucial 
to our mechanism. One might try to implement our idea 
replacing the vector-like doublets with SM singlets. This 
however would imply that the light neutrinos are partially 
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sterile which in turn would lead to non-universality in weak 
processes. For instance, the charged-current universality 
between /3-decay and /x-decay will be modified. As these 
observables are measured below the 0.1% level [|]], the lower 
bound on M is strengthen, significantly suppressing the BAU. 
Alternatively, one might study Type-II seesaw models in 
which Majorana neutrino masses are induced by the VEV of 
an electroweak triplet scalar boson; but constraints from the 
precision electroweak data are even more severe. 

Discussion. In this letter, we showed that the framework 
of late neutrino masses naturally leads to a viable model of 
electroweak baryogenesis that can be directly tested in near 
future experiments. 

The additional vector-like leptons with mass close to the 
EW scale play a crucial role in our scenario. This raises 
a coincidence problem since a vector-like mass term is 
unprotected by the symmetries of our model and therefore 
is unrelated to the EW scale, similar to the /i-problem in 
the Minimal Supersymmetric SM. There are however many 
proposed solutions to this problem: vector-like mass may 
be induced from supersymmetry breaking, additional Higgs 
fields, or strong dynamics. 

The first-order EW phase transition is mandatory. It may 
be induced, for instance, by higher order terms in the effective 
potential that arise from integrating out additional heavy 
states [Q. 

Our baryogenesis mechanism relies only on physics at 
the TeV scale, and hence is compatible with any scheme 
for new physics above the TeV scale. For example there 
could be flat extra dimensions just above the TeV scale, 
or a warped extra dimension as in the Randall-Sundrum 
scheme; in supersymmetric theories our mechanism allows 
for a low reheating temperature after inflation, solving the 
various moduli/gravitino problems. Finally, our scenario 
may accommodate unification by completing the vector-like 
leptons to full SU(5) representations. 
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